SOLUTION TO RUBEL'S QUESTION ABOUT DIFFERENTIALLY 
■ ALGEBRAIC DEPENDENCE ON INITIAL VALUES 
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Abstract. We prove that, for generic systems of polynomial differential equa- 
tions, the dependence of the solution on the initial conditions is not differen- 
tially algebraic. This answers, in the negative, a question posed by L.A. Rubel. 



■ 1. INTRODUCTION 

i-^h ' This work answers a question posed by L.A. Rubel in a paper which presents 

many research problems on differentially algebraic functions (see also ||). We recall 
that an analytic function f(z) (defined on some open subset of either R or C) is 
called Differentially Algebraic (DA) if it satisfies some differential equation of the 
form 

^| (1) Q(z,f(z)J'(z),...,fW(z)) = 0, 

^*T* | for all z in its domain, where Q is a nonzero polynomial of n + 2 variables. A 

f~L) ■ function of several variables is DA if it is DA in each variable separately when 

fixing the other variables. It is known that a function is DA precisely precisely 
CNJ ■ when it is computable by a general-purpose analog computer H . 

| Problem 31 of || asks: 

| Given a 'nice ' initial value problem for a system of algebraic equa- 

d • tions in the dependent variables yi, ■■■,y n , must yi(xo) be differen- 

tially algebraic as a function of the initial conditions, for each 

> 

X 



x ? 



■ In other words, let 



Rubel adds: "we won't say more about what 'nice' means except that the prob- 
lem should have a unique solution for each initial condition in a suitable open set" . 



( p ) Vk=Pk(yi,-,Vm), k = l...m 



be a system of differential equations, withp^ polynomials. Denote by yk{ri, ...,r m ;x) 
the solutions of (P) with initial conditions 

(2) y k (0) = r k , k = l...m. 

The question is whether it is true that the dependence of yi on rj (fixing x and 
the other initial conditions r k , k =/= j) is DA. 
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Here we first answer the question in the negative by constructing an initial value 
problem: 

(3) 2/i =2/2 - 2/1, 2/1 (o) = r x 

(4) 2/2=2/2(2/2-2/1), 2/2(0) = r 2 

and proving that, for most cc € M, both yi(0, r 2 ; a;) and 2/2(0, r 2 ; a;) are not DA with 
respect to r 2 . 

Moreover, we prove that this phenomenon is in fact generic in the class of non- 
linear polynomial differential systems (P) (and thus the answer to Rubel's question 
is very negative). To state this precisely, let S(m,d) denote the set of all poly- 
nomial systems (P) of size m, where the polynomials pi are of degree at most d. 
By identifying the coefficients of the various monomials occuring in the p^s with 
coordinates, S(m, d) is a finite dimensional vector-space, so it inherits the standard 
topology and measure. 

In all our statements the term generic will have the following meaning: when 
U C R N is an open set and V C U, we will say that V is generic in U if its 
complement in U is both of measure and of first Baire category. 

For each system P G S(m, d) and initial conditions r = (ri,...,r m ) G R m , let 
7(P, r) denote the largest interval containing for which the solutions of the initial- 
value problem (P),(|^) are defined. Let 

(5) A(P) = {(x,r) I r G W m , x G I(P,r)} 

(A(P) is the domain of definition of the functions 2/i(?"i, ■■■j?"m;2;) - we note that 
standard existence theory of ODE's implies that it is an open set; it is also easy to 
see that it is connected). For each 1 < j < m let 

^j( p ) = {fori, — ,rj-i,rj+i, ...,r m ) | 3rj such that (x,n, ...,r m ) £ A(P)}. 

Theorem 1. isswrae m > 2 and d > 2. For generic P £ S(m, d), we have: 

For any 1 < i,j < m, and for generic choice of (x, r%, Vj+x, r m ) G Aj(P), 

the function 

(6) fij(z) = yi(n, ...,rj-i, z,rj+i, ...,r m ;x) 
is not DA. 

We note that the construction of a specific example (||)-(f|) is essential for our 
proof of the 'generic' result of theorem [l]. 

The significance of the problem posed by Rubel is brought out when we note 
that the class of DA functions is a very "robust" one: 

(i) The class is closed under many of the constructions of analysis: algebraic oper- 
ations, composition of functions, inversion, differentiation and integration. 

(ii) The components of an analytic solution of a system of algebraic differential 
equations are themselves DA ||. 

(iii) A solution of a differential equation R(z, f(z), f'(z), f^ n \z)) = 0, where R 
is a DA function, is itself DA 0. We note a special case n = of this result, which 
can be called 'the DA implicit-function theorem', which we shall have occasion to 
use later: 

Lemma 2. Assume R(w, z) is a DA function, and f{z) is a real-analytic function 
satisfying 

R(f(z),z) = 
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for all z in some interval. Then f is DA on this interval. 

For the above reasons, most of the transcendental functions which are encoun- 
tered in "daily life" are DA (notable exceptions are the Gamma function and Ric- 
mann Zeta function, see ||). It is then natural to wonder whether functions ob- 
tained by the construction of looking at the dependence of the "final value" on 
the "initial value" (in other terms, the components of the "Poincare map" ) of an 
algebraic differential equation are also DA. Our results show that this is almost 
always not so. 

The following remarks show that the result of theorem |l| is in some sense close 
to optimal: 

(i) The condition d > 2 cannot be removed. Indeed if al — the solutions of the 
system are linear functions in x, and if d — 1 the solutions are linear combinations 
of exponentials in x, so the dependence on initial conditions is certainly DA. 

(ii) The condition m > 2 cannot be removed. Indeed if m = 1 then we are dealing 
with an equation of the form 

v =p(y), 

which is solved in terms of elementary functions, which implies that the dependence 
on the initial condition is DA. 

(iii) The restriction to generic systems cannot be removed, since, for example, when, 
for each k, depends only on y^, we are back to a decoupled system of equations 
of the same form as in (ii) above, so that again we have DA dependence on initial 
conditions. However, it might still be possible to obtain a stronger statement as to 
the size of the set of 'exceptional' systems. For example, we do not know whether 
the set of systems P 6 S(m, d) for which all the functions /y are DA is nowhere 
dense in S(m, d). 

In section || we prove a general result which underlies the proof of our 'genericity' 
results. In section || we construct some explicit examples of functions which we prove 
are not DA. In section || we show that the functions constructed in section || arise 
in the solution of the system ([|)-(|]), and this is used to prove that the solutions 
of this system are not DA in the initial conditions. Finally, in section ^| we use the 
specific example constructed in section ||, together with the general result of section 
|, to prove theorem [ij 



2. The Alternative Lemma 

In this section we present a result which we term the "alternative lemma" , which 
is used several times in the arguments of the following sections. This result says 
that if we have a parametrized family of analytic functions (with the dependence 
on the parameters also analytic), then either: (I) all the functions in the family 
are DA, or (II) generic functions in the family are not DA. Thus, to show that a 
generic function in a certain family is not DA, it is sufficient to prove that one of 
the functions in the family is not DA. 

The proof of the lemma is based on two facts: 
(i) The fact that an analytic function on a connected open set cannot vanish on a 
'large' subset without vanishing identically. 
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(ii) The Gourin-Ritt theorem (!]] which says that any DA function f(z) in fact 
satisfies a differential equation of the form ([!]) where Q is a polynomial with integer 
coefficients. 

We state the next lemma in both 'real' and 'complex' forms, since we shall have 
occasions to use both. 

Lemma 3. Let K = M or C. And let F : ft — > R, where il C W x K is open and 
connected, be an analytic function. For each u £ W , define 

fl(u) = {v e K | (u,v) e n} 

And let 

n' = {u g w I n(it) ^ }. 

TTien i/ie following alternative holds: either 

(I) For every u € fi', F(u,.) is DA on Q(u). In fact, all the functions F(u,.), 
u € f2', satisfy the same differential equation. 

or: 

(II) For generic u € Q', -F(it, .) is not DA on £l(u). 

PROOF: Let Q be the set of all polynomials with integer coefficients. Q is a 
countable set. For each Q G Q, let fl'(Q) be the set of u € f2', for which / = F(u, .) 
satisfies (0) on Q(u), We claim that either il'(Q) is of measure or Q'(Q) = ft' . 
Assume that il'(Q) has positive measure in R p . We choose open balls B\ C ft' and 
E>2 C K such that (i) B\ x B2 C and (ii) f2'(Q) n Si has positive measure. Then 
the real-analytic function 

g(u, v) = Q(v, F(u, v), D v F{u, v), D^F(u, v)) 

vanishes on the set (ft'(Q) PI Bi) x B2, which is of positive measure. By analyticity 
and the fact that Q is connected, this implies that g vanishes throughout Q, which 
means that alternative (I) holds. So if alternative (I) does not hold, f2'(Q) must be 
of measure 0. Since it is easy to see that ft'(Q) is a relatively closed set in 17', the 
fact that it has measure implies that it is nowhere dense. Since this is true for 
any Q £ Q, the countable union 

K = u QeQ n'(Q), 

is a set of measure and first category. By the Gourin-Ritt theorem, for all u 
outside K, F{u, .) is not DA. Thus (II) holds. 

3. Some functions which are not differentially algebraic 

In this section we define some new functions and prove that they are not DA. 
These results will be used in our construction of an explicit differential equation 
with non-DA dependence on initial conditions in section |ij 

We define the function H (c) (of c G C) by: 

(7) H(c)= ' < ! " 



For ceCto belong to the domain of definition Djj of H we need to ensure that 
the denominator ce u — u does not vanish for any u > 0. It is then easy to check 
that 

(8) D H = C-[0,-]. 

e 
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Lemma 4. The function H{c) is not DA. 

PROOF: We define h(z) — H(^). We shall show that h is not DA, which implies 
that H is not DA. We expand h in a power series (which converges for \z\ < e): 

/ ze du / ze _ u y z k u k e -ku du 

In 1 - ZUe U In 



fc=0 
oo 



fe=0 ^° fc=l 

We now use the theorem of Sibuya and Sperber pj, which gives the following 
necessary condition on a power series X^fcLo a fc 2 with rational coefficients in order 
for it to satisfy an algebraic differential equation: 

(9) \a k \ P <e ck forallk>0, 

where \a\ p (p prime) is the p-adic valuation defined by writing a — p 1 ^ (m, n not 
divisible by p) and setting \a\ p = p~ l . Here we use this necessary condition with 

p = 2. We have at — ^pr^-- We now choose k = 2 J . The heighest power of 2 
dividing (k — 1 )! is less than 

fc-lfc-lfc-1 

— + — + — + ...<k-l = »-l. 
On the other hand the denominator of a& is 2 j2J . These two facts imply that 

which together with (^) implies: 

2^- 1)23+1 <\a 23 \2<e C2 \ 

or, taking logs 

((i - 1)2* + l)/o fl (2) < CT, 

which is obviously false for j sufficiently large. Hence the necessary condition of 
Sibuya and Sperber is not satisfied for our series, completing the proof. 

We note a few properties of the restriction of H to the real line (which is defined 
outside the interval [0, -]), obtained by elementary arguments. See figure |l| for the 
graph of H , plotted using MAPLE. 

(i) H is decreasing on oo), with lim c ^i + H{c) — +oo, lim^oo H{c) = 0. 

(ii) H is decreasing on (— oo, 0), with lim c _ ( _ 00 H(c) = 0, lim c ^ - H{c) = — oo. 

We now define the function F(s, c), where s 6 1 is considered a parameter and 
c e C by: 

du 



(10) F(s,c)-- 

Since for F(s, c) to be defined we need to ensure that the denominator does not 
vanish for u £ [0, s], it can be checked that the domain of definition of F is 

Dp = {(s, c) | s £ M, c £ C - J(s)}, 
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Figure 1. Graph of the function H(c) defined by (f7j) 



where 

J(s) = [se" s ,0] for s < 
J(s) = [0, se- s ] for < s < 1 

J(s) = [0, -] /or s > 1. 

e 

Using the fact that H is noi DA (lemma ||), we now prove that the function 
F(s, .) is not DA for generic s. The argument is similar to the one in 0] showing 
that the function z — > u z e~ u du is not DA for generic a, 6 by using the fact that 
the Gamma function is not DA. 

Lemma 5. For generic set, the function F(s, .) is not DA. 

PROOF: We apply lemma ||| with p = 1, K = C, fi = Dj? clxC (which is indeed 
open and connected). We have ft' = M. We want to show that alternative (II) of 
lemma || holds. Let us assume by way of contradiction that (I) holds, so that F(s, .) 
is DA for all s£l, and in fact there is a common polynomial differential equation 
satisfied by all F(s, .): 

(11) Q(c,F(s,c),D c F(s,c),...,D^F(s,c)) = 0, for all (s,c) e D F . 
We now note that 

lim F(s, c) = H(c) 

s — >oo 

uniformly on compact subsets of Dh- This implies that derivatives of D^F(s,c) 
also converge to corresponding derivatives of H, and thus from ([ll]) it follows that 

Q(c,H(c),H'(c),...,H^(c))=0 
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for all c e Dh, hence H is DA, in contradiction with lemma |]. This contradiction 
proves that alternative (II) holds, as we wanted to show. 

We remark that we would guess that the 'exceptional' set in lemma || is {0}, so 
that in fact F(s, .) is not DA for all s =/= 0, but we do not know how to prove this. 

We now restrict c to be real nonzero number, and we define a new function 
G(t, c) by the relation: 

(12) G(F(s, c), c) = s for all (s, c) G D F . 

In other words we now look at c as a parameter and define G(.,c) as the inverse 
function of F(.,c). To see that G is well-defined and determine its domain of 
definition, we note the following properties of F, which are elementary to verify. 

Lemma 6. For each c < - let s*(c) denote the solution (in the case < c < -, 
the smaller solution) of the equation se~ s — c. We have: 

(i) When c < the function F(.,c) is decreasing on (s*(c), oo), 
lim s ^ s »( c ) + F(s, c) = +oo and lim^oo F(s, c) = H(c) < 0. 
Hence G(.,c) is defined on (H(c), +oo). 

(ii) When < c < - the function F(.,c) is increasing on (— oo, s*(c)), 
linig^-oo F(s, c) = — oo and lim s ^ s *( c )_ F(s, c) = +oo. 

Hence G(.,c) is defined on (—00,00). 

(Hi) When c > i the function F(., c) is increasing on (—00, 00), 
lim^-oo F(s, c) = —00 and lim^oo F(s, c) = H (c) > 0. 
Hence G(.,c) is defined on (— 00, H(c)). 

Putting these facts together, we obtain that the domain of definition Dq of G is: 

D G = D G UD+, 

where 

D G = {(t, c) I c < 0, H(c) < t < 00}, 

D± = {(t,c) I < c < -} U {(t,c) I c> - -00 < t < H(c)\. 

e e 

( in figure Dq is the domain bounded by the left part of the graph of H and the 
t-axis, and Dq is the domain bounded by the right part of the graph and the t-axis). 

We also note a fact that seems hard to prove directly, but which follows indirectly 
from the results of the next section, as will be pointed out. 

Lemma 7. The function G can be continued as a real-analytic function to the open 
connected domain 

(13) D = D G UD+ U {(i,0) I t e M} 
by setting 

(14) G(i,0)=0 forallteR 

(in figure^ D is the domain bounded between the two parts of the curve representing 
the graph of H). 

Lemma 8. For generic f G R, the function G(t, .) is not DA. 
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PROOF: Assume by way of contradiction that the conclusion of the lemma is not 
true, so that G(t, .) is DA for a set of values of t which is either of positive measure 
or not of first category. Then, applying lemma ||, with Q = D (recall (|l3|)), we 
conclude that G(t, .) is DA for all t € R. We note also that G(.,c) is DA for all 
value of c, since for c 7^ it is the inverse of F(., c) which is DA, and the inverse 
of a DA function is DA, while for c = it is identically by Thus G is DA 

with respect to both variables on Q,. Now fixing any s € R, and using the defining 
relation ( jl^ ) we have 

G(F(s, c),c) - 3 = 

for all c. We now use lemma ||, with R(w, z) — G(w, z) — s, to conclude that F(s, .) 
is DA. Since s £ R was aribtrary, we have that F(s, .) is DA for all s e R, but this 
contradicts the result of lemma |^. This contradiction concludes the proof. 

4. Construction of a differential equation with non-DA dependence 

on initial conditions 

Theorem 9. Consider the initial value problem: 

(15) y[ = y 2 - yi(0) = n 

(16) 2/2 = V2(y2 - 2/1), 2/2(0) = r 2 

Let the solutions be denoted by yi{r\,r2]x), 2/2 (t, ^2; %), defined on the maximal 
interval I(rt, ^2)- Then for generic i£K, the functions 2/1 (0, r^x), 2/2(0, r 2 ; x) are 
not DA with respect to r 2 . 

To prove theorem ^, we assume for the moment that r 2 ^ 0, which implies that 
2/2(2;) ^ for all x (since by ( |l6| ) and by the standard uniqueness theorem for 



ODE's, if 2/2 vanishes at one point, it vanishes identically), and note that from (15) 
and ( |l6| ) we have 

2/2 

which implies 

V2{x) =y 2 (o) e ^-y^\ 

and in the case 2/1 (0) = r\ = 0, 2/2(0) = r 2 : 

(17) y 2 (x)=r 2 e y ^ x \ 

Substituting (^) back into (|l5|), we have 

yi(ar) = r 2 e yi(;r) - j/i(x), 

or, 

r2e yi(a;) _ y 1 (a;) 

Thus, integrating, and recalling the function F(s,c) defined by (|l0|), we have: 

F{yi{0,r 2 ;x),r 2 ) = x 
for all x € 7(0, r 2 ). In other words, recalling the definition (O) of G(f , s), we have: 



(18) 



2/i(0, r 2 ; x) = G(x, r 2 ) /or r 2 /0,i(E 7(0, r 2 ) 
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We note in passing that, since yi(0, r2]x) is real-analytic in its domain of defi- 
nition, and is defined also for T2 — 0, with yi(0, 0;a;) = 0, ( |l8| ) implies that G can 
be extended to D as lemma [?] claimed. 

In lemma we showed that G(t, c) is not DA with respect to c for generic t. 



From (18) we then get, for generic x, that y\{Q 1 r2',x) is not DA with resepct to r<2- 
Using this fact and (|l7|), together with the fact that compositions and products of 
DA function are DA, we get the same conclusion for , concluding the proof of 
theorem |9| 

5. PROOF OF THE MAIN THEOREM 

We first set some notation. Define E C S(m, d) x M x R m by 

Y, = {{P,x,r) | PeS(m,d), rel™, x€l(P,r)}. 
For each 1 < j < m let 

E^- = {(P, x, ri, rj-x, r 3 -+i, r m ) | 3rj suc/i t/iai (P, cc, r) e E}. 
The proof of theorem [|, will be based on the following 
Lemma 10. Assume m > 2, d > 2. Fix 1 < j < m. 

Then for generic (P, n, fy-i, Tj+i, ...r TO ,:c) € E^, the functions (1 < i < 
defined by (f^jj are not £L4. 

PROOF: We note first that we may, without loss of generality and for notational 
convenience, assume that j = m. Define Fk : E — > M (1 < A; < m) by 

F fe (P,x,ri, ...,r m ) = y/s(ri, ...,r m ;a;). 

These are real analytic-function on E, and we may apply lemma |^, with K = K, 
fi = E, u identified with (P, x, r\, ...,r m _i), v identified with r m , and F being any 
of the Fk (1 < fc < m). We would like to rule out alternative (I). We now fix some 
1 < i < m — 1, and note that among the systems P 6 S(m, d) is the system 

J/i = Vrn - Vi 

y'm = ym(y m - y») 

y*: = for i,m. 

By theorem |^, for this system, for = and for arbitrary values of (k ^ i, m), 
the dependence of yi(rx, r m ; x) and of y m (ri, r m ; x) on r m is not DA. This 
rules out alternative (I) of lemma || both for P, and for F m , which implies that 
alternative (II) holds for both. Since 1 < i < m — 1 is arbitrary, this means that 
for each 1 < k < m there is a generic subset E| of Yj' m for which the function Fk is 
not DA with respect to r m . Thus, setting 

e 9 = n x < fe < m ES 

we get a generic set, and for each (P, x,r\, ...r m -x) G E 9 we have the conclusion of 
lemma ^0|. 

Finally, to derive theorem |l| from lemma [l0|, we recall two classical results (see 
H). Let C C R N x K M , and define for each u G R N , 

C(u) = {u G M M | («,«) G C}. 

We have: 
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(i) If the set C is of measure then C(u) is of measure in M M for all u £ M. N 
except a set of measure 0. This is (a special case of) Fubini's theorem. 

(ii) If the set C is of first category then C(u) is of first category in R M for all 
u G except a set of first category. This is the Banach-Kuratowski theorem. 

PROOF OF theorem [l]: We will prove the conclusion of theorem |l| for all P £ 
S(m, d) outside a set of measure 0, making use of Fubini's theorem. To prove the 
same conclusion for P outside a set of first category, one only needs to replace 
Fubini's theorem with the Banach-Kuratowski theorem. 

For 1 < i,j < m define C,j to be the set of (P,x,ri, Tj+i, ■■■,r m ) € E^- 
for which the function defined by @ is DA. By lemma |l^, each of the sets Cy is 
of measure and of first category. Hence by Fubini's theorem, for all P € S(m, d) 
except a set of measure 0, the set 

c %j( p ) = {( x ,n, -,rj-i,rj + i, ...,r m ) \ {P,x,r x , r,_i, r 3+ i, ...,r m ) e C lj } 

has measure 0. Below we show that the sets Cij(P) are also of first category for 
any P for which they are of measure 0. Since this conclusion is valid for any 
1 < i,j < m , we have that for all P £ S(m, d) except a set of measure 0, all the 
sets Cij (P) have measure and are of first category. This is precisely the content 
of theorem [l[ 

It is left to show that whenever dj(P) is of measure 0, it is also of first category. 
Without loss of generality and for notational convenience, we assume j = m. We 
apply lemma H with K = R, fi = A(P) (recall (||)), u identified with (x,ri, ...,r m _i), 
v identified with r m , and F — fi m . Since we assume that Ci m (P) is of measure 0, 
so that F is not DA with respect to r m for almost all (x, n, r m _i) G A m (P), 
alternative (I) of lemma || certainly does not hold, so (II) holds which means that 
Cim(P) is indeed also of first category, concluding the proof. 
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